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1.  SUMMARY 


It  is  known  that  friction  due  to  boundary  layer  disturbances  can  affect  the  performances  of 
many  engineering  systems  and  there  are  many  works  on  the  suppression  of  the  transitional 
boundary  layer.  Most  of  the  results  in  the  literature  are  for  the  flat-plate  boundary  layer 
but  the  behaviour  of  the  velocity  and  pressure  changes  with  the  curvature.  This  work  aims 
to  extend  the  results  of  the  flat-plate  boundary  layer  to  a  Rankine-body  boundary  layer. 

A  small  amplitude,  unsteady,  three-dimensional  perturbation  of  a  mean  base  flow  is 
considered.  The  flow  is  taken  incompressible  and  the  Reynolds  number  high.  The  analysis 
is  performed  in  a  new  special  coordinate  system,  where  the  streamwise  coordinate  is  replaced 
by  the  potential  function  and  the  wall  normal  coordinate  by  the  stream  function. 

The  external  flow  is  decomposed  into  a  mean  flow  and  a  perturbation  flow.  The  first 
one  is  known  and  it  is  given  by  a  uniform  flow  which  encounters  a  source  flow  because  the 
presence  of  the  body  can  be  seen  as  a  point  source.  The  perturbation  flow  around  and 
far  from  the  body  is  solved  by  applying  WKBJ  theory,  a  method  to  obtain  approximated 
solutions  of  linear  differential  equations  whose  the  highest  derivative  is  multiplied  by  a  small 
parameter. 

In  the  viscous  region  near  the  body,  the  flow  is  also  decomposed  into  a  mean  flow  and  a 
perturbation  flow.  The  base  flow  satisfies  the  steady  mean  flow  equations  in  optimal  coor¬ 
dinates  and  since  it  depends  on  the  streamwise  coordinate  there  is  no  similarity.  Boundary 
region  equations  in  optimal  coordinates  for  the  perturbation  flow  are  obtained.  The  bound¬ 
ary  region  equations  are  numerically  solved  by  imposing  no-slip  conditions  at  the  wall  for 
the  three  velocity  components  and  appropriate  boundary  conditions  at  a  large  distance  from 
the  body.  The  latter  are  found  by  asymptotic  matching  between  the  inner  and  the  outer 
solution.  Differently  from  the  flat  plate  geometry,  the  outer  pressure  disturbance  has  a  role 
on  the  boundary  layer  dynamics  and  the  inner  perturbation  pressure  matches  the  outer  one 
at  a  large  distance  from  the  body.  Another  important  difference  from  the  flat  plate  case 
is  that  there  is  a  pressure  gradient  in  the  inner  mean  flow.  Boundary  region  equations  are 
parabolic  in  the  streamwise  coordinate  thus  a  marching  procedure  along  this  coordinate 
will  be  used  in  the  code.  Initial  conditions  to  start  the  code  are  derived. 

The  results  suggest  that  the  mean  pressure  gradient  modifies  the  shape  of  the  streamwise 
velocity  profile  compared  to  the  one  of  the  flat-plate  boundary  layer.  It  is  observed  that 
the  streamwise  wavenumber  plays  a  key  role  in  the  streaks  development. 
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2.  INTRODUCTION 


The  interaction  of  a  fluid  in  motion  with  a  solid  surface  affects  the  performance  of  mechan¬ 
ical,  aeronautical,  civil  and  chemical  engineering  systems.  When  an  object  moves  through 
a  fluid,  a  boundary  layer  is  created  on  its  surface.  In  this  region  viscous  effects  are  not 
negligible.  Boundary  layer  can  either  be  laminar  or  turbulent.  In  the  first  one  the  fluid 
flows  in  parallel  and  ordered  layers  while  in  the  second  one  the  flow  is  chaotic.  Generally, 
the  boundary  layer  changes  from  laminar  to  turbulent  through  a  transitional  zone.  It  has 
been  established  that  turbulent  flows  exert  a  much  higher  wall  friction  than  laminar  flows. 
Since  vehicles  such  as  cars  and  airplanes,  for  example,  consume  an  enormous  amount  of 
energy  due  to  friction,  many  works  have  been  directed  to  the  suppression  of  transitional 
boundary  layer  disturbances.  It  is  better  achieving  drag  reduction  in  the  still  laminar 
transitional  zone  rather  than  in  the  already  developed  turbulent  one  because  suppressing 
small-amplitude  laminar  disturbances  requires  less  energy.  The  reduction  of  drag  has  a  great 
technological  importance  because  it  implies  a  significant  decrease  of  the  enormous  amount 
of  energy  consumed  by  airplanes  during  flight,  moreover  flight  costs  and  aerodynamic  noise 
could  be  reduced  and  number  of  passenger  per  flight  could  be  increased.  Such  outcomes 
could  limit  the  global  warming  effect.  Despite  the  process  which  governs  the  transition  has 
been  widely  studied,  at  present  is  not  fully  understood. 

The  flow  within  the  boundary  layer  strongly  depends  on  external  flow.  It  is  therefore  of 
utmost  importance  to  study  the  effects  of  free  stream  turbulence  (FST)  on  the  boundary 
layer.  FST  are  defined  by  the  turbulence  level  {Tu),  the  turbulence  length  scales  and 
the  turbulence  spectrum.  Experimental  works  have  confirmed  that  if  Tu  is  low  (  <  1  % 
)  the  transition  from  laminar  to  turbulent  boundary  layer  occurs  due  to  the  growth  of 
Tollmien-Schilichting  (TS)  waves;  otherwise,  if  Tu  is  high  (  >  1%)  the  mechanism  which 
determines  the  transition  from  laminar  to  turbulent  is  called  bypass  transition.  During 
bypass  transition  the  flow  is  laminar  and  it  is  characterised  by  long  streamwise  perturbations. 
These  streaks  are  also  called  Klebanoff  distortions,  because  Klebanoff  was  among  the  first 
to  observe  them.  Due  to  non-parallel  flow  effects,  these  low-frequency  components  of  three- 
dimensional  vortical  disturbances  in  the  free  stream  penetrate  the  boundary  layer  and  they 
produce  significant  distortions  in  the  spanwise  direction.  There  is  a  large  number  of  direct 
laboratory  investigations  of  the  transition  process  in  the  presence  of  Klebanoff  fluctuations. 

?  observed  intermittent  appearance  of  wave  packet  in  the  boundary  layer  at  moderate 
levels  of  free-stream  turbulence.  Even  if  he  did  not  understand  the  origin  of  Klebanoff 
modes  he  performed  many  experiments  which  revealed  some  important  characteristic  of 
the  wave  packets.  They  appear  only  when  the  turbulence  level  is  higher  than  a  threshold, 
they  grow  faster  than  TS  waves  and  they  propagate  quickly  downstream  and  slowly  in 
the  spanwise  direction.  Downstream  the  streaky  region  turbulent  spots  appear  and  then 
turbulent  boundary  layer  fully  develops.  While  TS  wave  transitions  represent  idealised 
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situations,  bypass  transitions  describe  real  situations  for  industrial  flow  systems. 

Bypass  transition  has  been  largely  studied  experimentally  and  all  the  results  showed 
that  FST  induced  streamwise  streaky  structures  of  alternate  high  and  low  velocity  in  the 
boundary  layer.  ?  carried  out  experiments  using  techniques  like  flow  visualisation  and 
hot-wire  anemometry  in  order  to  determine  the  effects  of  FST  (in  the  range  of  1-6  %)  in  the 
boundary  layer.  They  performed  the  experiments  using  grids  of  various  sizes  to  generate 
the  free-stream  turbulence  on  a  test  plate.  They  found  that  the  spanwise  scale  decreased 
moving  downstream,  where  it  approached  the  boundary  layer  thickness.  Conversely  the 
streamwise  scale  increased  and  it  became  proportional  to  the  boundary  layer  thickness, 
going  downstream.  They  also  pointed  out  that  the  breakdown  of  the  streaky  structures  was 
associated  with  a  secondary  instability  which  seemed  to  lead  to  turbulence. 

Many  authors  performed  direct  numerical  simulations  (DNS)  on  bypass  transition  in¬ 
duced  by  free-stream  turbulence  (sometimes  in  complement  with  laboratory  experiments).  ? 
depicted  eddies  and  interactions  during  transition  induced  by  free-stream  turbulence.  Their 
simulations  showed  that  when  low-frequency  modes  penetrated  the  boundary  layer,  they 
produced  lower  frequency  modes,  which  were  amplified  by  the  shear  and  they  developed 
into  streaks  in  the  streamwise  velocity.  Streaks  were  the  same  of  the  experimental  measure¬ 
ments  of  Klebanoff.  Simulations  of  ?  were  based  on  Orr-Sommerfeld  theory.  This  choice 
allowed  to  use  a  rectangular  domain  and  since  the  lower  wall  extended  over  all  domain,  a 
better  algorithm  could  be  used.  Results  showed  that  streaks  were  stable  if  they  remained 
near  the  wall;  instead  if  they  lifted  up,  the  disturbance  rapidly  grew  and  turbulent  spots 
appeared.  Then  the  spots  led  to  turbulent  boundary  layer.  ?  explored  the  penetration 
of  external  vortical  disturbances  in  the  boundary  layer,  by  studying  Orr-Sommerfeld  and 
Squire  modes  interaction.  The  evolution  of  modes  was  analysed  by  DNS  and  the  domain 
included  boundary  region  zone.  The  transition  process  was  simulated  with  two  modes:  one 
strongly  coupled  and  one  weakly  coupled  high-frequency  mode.  Results  showed  that  if  one 
mode  penetrated  and  the  other  one  no,  there  was  transition.  In  this  paper  it  was  confirmed 
that  low  frequencies  penetrated  the  boundary  layer  and  they  led  to  Klebanoff  modes. 

For  a  long  time  bypass  transition  has  only  been  experimentally  observed,  but  it  could 
not  be  explained  theoretically  .  The  optimal-growth  theory  was  first  developed  by  ?.  The 
objective  of  the  analysis  is  to  find  the  optimal  perturbations,  namely  the  initial  velocity 
profile  which  produces  the  biggest  amplification  of  the  disturb.  In  the  Orr-Sommerfeld 
theory  the  continuous  spectra  of  the  Orr-Sommerfeld  equations  is  used  to  describe  the 
streaks  generated  by  free-stream  turbulence. 

?  analysed  the  effects  of  free-stream  turbulence  on  the  pretransitional  boundary  layer 
in  a  flat  plate.  The  free-stream  flow  is  decomposed  in  a  mean  flow  and  a  purely  convected 
three-dimensional  perturbation  flow.  Their  results  showed  that  the  free-stream  streamwise 
vorticity  is  responsible  for  the  generation  of  the  streamwise  streaks  inside  the  boundary 
layer.  The  streaks  are  initially  amplified  then  they  decay  downstream. 

?  generalised  the  results  of  ?  to  the  unsteady  case.  They  studied  the  instability  of 
unsteady  nonlinear  streaks  due  to  free-stream  disturbances  on  an  incompressible  laminar 
boundary  layer.  Their  attention  was  focused  on  long  disturbance  wavelengths  because  it 
is  known  that  they  are  responsible  for  the  streaks  generation.  Their  results  showed  that 
the  nonlinearity  can  inhibit  the  amplification  of  streaks,  thus  it  has  a  stabilising  effect  and 
nonlinear  interactions  have  an  effect  of  distortion  on  the  mean  flow.  These  numerical  results 
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are  in  agreement  with  the  experimental  results  of  ?. 

?  focused  on  the  solution  of  the  external  flow  around  a  body  of  generic  shape  where 
a  small-amplitude,  unsteady,  viscous  motion  was  imposed  on  steady  potential  flow.  He 
derived  the  equations  for  the  perturbation  velocity  and  then  divided  the  solution  for  it 
in  two  parts:  the  first  one  is  the  homogeneous  solution  which  satisfies  the  momentum 
equation  without  the  forcing  term  (namely  without  the  pressure  gradient);  the  second  one 
is  the  potential  solution,  from  which  the  outer  pressure  is  derived.  The  total  perturbation 
velocity  is  given  by  the  summation  of  these  two  solutions.  He  first  found  the  homogeneous 
solution  by  WKBJ  theory  and  then  from  the  continuity  equation  he  calculated  the  potential 
solution.  By  using  this  theory  he  looked  for  the  phase  and  the  amplitude  of  the  perturbation 
velocity.  The  equation  for  the  phase  in  a  new  coordinate  system  is  found  to  be  the  Eikonal 
equation.  He  found  a  complete  analytic  solution  of  the  perturbation  velocity  by  changing 
coordinates.  ?  supposed  that  the  steady  potential  flow  is  not  influenced  by  viscous  effects 
and  to  study  the  problem  he  included  viscous  effects  in  rapid  distortion  theory.  He  used 
his  results  to  calculate  the  flow  over  a  circular  cylinder  and  he  showed  the  numerical  results 
for  large  Reynolds  numbers.  In  particular  he  focused  on  the  behaviour  of  the  streamwise 
perturbation  component  velocity  along  the  mean  flow  stagnation  streamline.  Its  absolute 
value  vanishes  when  the  dimensionless  wavenumbers  along  the  perpendicular  directions 
go  either  to  zero  or  infinity.  Moreover,  plots  of  the  same  function  versus  the  streamwise 
coordinate  show  that  it  decreases  as  the  coordinate  goes  to  zero. 

?  worked  on  rapid  distortion  theory,  to  extend  the  work  of  ?.  He  considered  compressible 
flows,  with  no  restrictions  on  the  Mach  number  value  with  high  Reynolds  numbers  and  he 
started  from  the  linear  inhomogeneous  wave  equation  of  a  potential  function  associated 
with  the  perturbation  velocity.  He  assumed  that  the  upstream  disturbance  was  a  small 
perturbation  of  the  mean  flow.  He  showed  that  the  perturbation  velocity  consisted  of  two 
parts:  one  was  a  known  function  of  the  imposed  upstream  distortion  and  the  mean  flow,  it 
had  zero-divergence  and  it  was  decoupled  from  the  fluctuations  in  pressure;  the  other  one 
was  an  irrotational  disturbance,  given  by  the  gradient  of  the  potential  function  which  was 
the  solution  of  the  wave  equation  and  it  was  related  to  pressure  fluctuations.  In  the  paper 
is  also  remarked  that  even  without  distortion  of  the  mean  flow,  pressure  fluctuations  could 
also  be  produced  by  boundary  conditions  on  the  surface  of  the  obstacle. 

The  effect  of  the  curvature  combined  with  a  favourable  pressure  gradient  on  the  streaks 
are  unknown. 


2.1  Objectives 

•  Computing  the  outer  flow  dynamics,  which  is  crucial  to  derive  boundary  conditions 
for  boundary  region  equations. 

•  Analyzing  the  effect  of  pressure  gradient  on  the  pre-transitional  Rankine  body  bound¬ 
ary  layer. 

•  Studying  the  evolution  of  Klebanoff  modes  inside  the  boundary  layer  and  identifying 
the  streaks  with  accelerating  mean  flow. 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


3.  METHODS,  ASSUMPTIONS  AND  PROCEDURES 


3.1  Mathematical  formulation  of  the  flow  around  a  Rankine  body 

The  flow  around  a  Rankine  body  is  considered.  A  sketch  of  the  geometry  is  shown  in  figure 
3.1.  A  small  amplitude,  unsteady,  three-dimensional  perturbation  flow  is  superimposed  to 
a  uniform  mean  flow  upstream.  The  body  can  be  represented  by  a  flow  produced  by  a 
point  source.  When  the  flow  encounters  the  body  both  the  mean  and  the  perturbation  flow 
are  distorted.  The  flow  returns  parallel  downstream,  it  is  assumed  incompressible  and  high 
Reynolds  numbers  based  on  the  characteristic  length  of  the  body  L*  are  considered  .  The 
domain  is  divided  in  two  main  asymptotic  regions:  the  outer  region  around  the  body  but 
distant  from  it  and  the  inner  small  viscous  boundary  layer  region  alongside  the  body.  The 
following  hypothesis  are  made: 

•  A*  <C  L*  where  A*  represents  all  the  three  wavelengths  of  the  disturbance.  This 
hypothesis  is  valid  upstream  and  it  will  be  used  during  the  analysis  of  the  outer 
perturbation  flow. 

•  The  mean  flow  is  steady  and  two-dimensional  therefore  periodicity  is  assumed  along 
z  and  t. 

•  The  solution  is  linear,  namely  e / k<[>  <C  1  inside  the  boundary  layer,  therefore  terms  of 
order  e2  can  be  neglected  in  the  equations. 

Table  3.1  gives  an  overview  of  the  asymptotic  regions  around  the  Rankine  body.  In  the 
following  sections  flows  in  the  outer  and  inner  region  will  be  analysed. 

3.1.1  Scaling  and  change  of  coordinates 

In  this  section  the  appropriate  scaling  is  developed  and  the  change  of  coordinate  is  described. 
Lengths  which  are  scaled  by  L*  (the  characteristic  length  of  the  body)  are  indicated  by  the 
subscript  L.  Lengths  which  are  scaled  by  A*  (the  spanwise  disturbance  wavelenght)  do  not 
have  any  subscript.  The  velocity  components  are  made  dimensionless  by  ,  pressure  by 
p*U^  where  p*  is  the  constant  density  and  time  by  L* /U^.  All  the  dimensional  quantities 
are  marked  by  the  symbol  *. 

The  following  continuity  and  Navier-Stokes  equations  have  to  be  satisfied  by  the  mean 
flow 


V-  V  =  0, 


(3.1) 


<9V  1 

+  V  •  VV  =  -Vp  +  — V  •  (VV), 
<76  L  -6tL 


(3.2) 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


3  ^,+8 

n  ^ 
c  n,> 


g.  3.1:  Perturbed  flow  around  the  body.  The  origin  of  Cartesian  and  polar  coordinate  sys 
is  in  the  point  source.  At  the  stagnation  point  a  rotate  coordinate  system  is  used, 
asymptotic  regions  around  the  body  are  shown. 
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Mean  flow 

Perturbation  flow 

Outer 

•  Inviscid-2D. 

•  Uniform  plus  source  flow. 

•  Pressure  gradient. 

•  Stagnation  point:  Hiemenz 
flow 

•  Inviscid-3D. 

•  WKBJ  method. 

•  Eikonal  equation. 

•  Equation  for  the  amplitude  velocity. 

•  Pressure  gradient. 

Inner 

•  Viscous-2D. 

•  Steady  nonlinear  mean  flow 
equation  in  optimal  coordi¬ 
nates. 

•  Viscous-3D. 

•  Boundary  layer  equations. 

•  Asymptotic  matching  with  outer 
flow  to  obtain  boundary  conditions. 

•  Derivation  of  initial  conditions. 

Tab.  3.1:  Decomposition  of  the  flow  around  the  Rankine  body 


where  i?L  =  U^L* /v*  3>  1  is  Reynolds  number  based  on  L* ,  the  V  operator  is  in 
coordinates  (xL,yL,zL)  and  V  is  the  vector  which  represents  the  total  mean  flow  velocity. 

Since  the  body  is  not  a  flat  plate,  it  is  convenient  to  work  in  a  special  coordinate  system 
where  the  streamwise  coordinate  ,xL  is  replaced  by  the  outer  mean-flow  velocity  potential  <1>L 
and  the  wall  normal  coordinate  yL  is  replaced  by  the  outer  mean- flow  stream  function  \I/L. 
In  this  new  coordinate  system  the  geometry  is  simplified  because  the  streamlines  follow  the 
shape  of  the  body,  therefore  the  body  is  a  flat  plate  in  the  new  system.  These  coordinates 
are  called  optimal  coordinates  and  in  Appendix  A  continuity  and  Navier-Stokes  equations 
in  optimal  coordinates  are  derived.  They  are 


~  Sicvt,  -dwq,  8wz 

U — -  -I-  U — -  -I - -  =  0 

STL  +  STL  +  dzL 


(3.3) 


1  dp  1 

TP  S $L  +  rl 


1  dw$  dw$  dw$  1 

U2  SfL  S«hL  STL  u  dzL 


dw $  Slog  U  ,  2  2  \  _ 


S^T 


d2w<$,  i  d2w<s>  i  1  d2w<$,  ^Sln  U  ( dwq,  dw$\  2  Sin  U  dwz 


d<5>2 


+ 


ST^  jj2  dz'l 


STl  V  Sd*L  StJ  u  S<hL  SzL 


(3.4) 
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l  dp  1 
U2  d^>L  +  Rh 


1  dwq,  dwq, 

+  'u;<i>  vrv — I-  '! 


d2W\ $ 

cWf 


U 2  dth 

d2w q. 


+ 


dT2 


+ 


d<S>  L 

1  d2iuq, 

U 2  dz 2 


+  2 


dwq/  1  dvjxji  d  log  U  ,  o  2  \ 
<9TL  u  dzh  <9TL  v  w 
d  In  U  (  dwq,  dwq,  \  2  d  In  U  dw- 

<9<hL  V  d$L  d^L )  u  d^h  dzL 

(3.5) 


1  dp  1 

U 3  dz^  Rl 


d2 


w- 


d$l 


1  dwz 

u2  dtL 

d2uh 

+ 


<9T2 


+  wq> 


+ 


dwz 

d$L 

1  d2wz 
U 2  dz 2 


drvz  1  dwz 

+  Wq,  — -  +  ^wz— — 

u  dzL 


wz 

+ 

17 


92[/  d2C/ 
+  5^2 


+ 


+  2 


9  log  [7  d  log  U 
- - h 


V 

dwz  d  In  U 
~d¥  <9<h T 


+ 


/ 

dwz  dlnU 
d%_lRK  ^ 
(3.6) 


where  U  =  |U|  is  the  amplitude  of  the  outer  mean  velocity  which  has  its  only  component 
along  the  streamlines,  w$  =  uq,/U,  wq,  =  uq, fU  and  wz  =  w/U,  with  and  uq,  the  velocity 
components  in  the  $  and  T  directions  respectively. 


3.1.2  The  mean  Sow 
Outer  mean  Sow 

As  explained  in  ?  a  combination  of  uniform  flow  with  a  source  is  well  represented  in  a  polar 
coordinate  system. 


r  =  \J  x2  +  y2, 


9  =  at  an 


t  ■  a  M 

VT  =  r  sm  9  -\ - 

2vr 


A 

2’ 


(3.7) 

(3.8) 


r  cos  9  + 


—  In  r  — 
2n 


A 
27 r 


(3.9) 


where  d>L  =  &* / (U^L*)  is  the  potential  function,  \HL  =  ^*/(U^0L*)  is  the  stream 
function  and  A  is  the  dimensionless  strength  of  the  source. 

The  absolute  value  of  the  mean  outer  velocity  can  be  analytically  found  as  a  function 
of  r  and  9. 


U  = 


A  cos  9 

1  + - + 


A2 


(3.10) 


nr  47t2?’2  ’ 

and  therefore  the  potential,  stream  function  and  the  mean  velocity  are  known  as  func¬ 
tions  of  r  and  9.  It  is  not  possible  to  analytically  invert  the  relationships  between  r-9  and 
T-T,  thus  this  is  numerically  performed  by  using  Newton-Raphson  algorithm. 

The  expression  of  r  as  function  of  9  and  TL  is 
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1 

r  =  - 

sin# 


A6  A 

TL - f-  — 

2ir  2 


The  asymptotic  behavior  of  the  inner  mean  flow  near  the  stagnation  point  in  optimal 
coordinates  along  the  body  (streamline  equal  to  zero)  has  been  found  by  ?  and  it  is 


uH  — 


(3.11) 


Figure  3.2  shows  that  in  the  neighbourhood  of  the  stagnation  point  the  asymptotic 
solution  and  the  mean  velocity  overlap.  Figure  3.3  shows  the  outer  mean  velocity  on  the 
body  surface  14  as  a  function  of  the  coordinate  s,  which  is  the  coordinate  along  the  body, 
defined  as 


s  = 


dye 

dxL 


dxL. 


Fig.  3.2:  Comparison  between  the  outer  mean  ve¬ 
locity  on  the  body  and  the  asymptotic 
solution  close  to  the  stagnation  point. 
In  the  neighbourhood  of  the  stagnation 
point  the  behavior  is  the  same. 


Fig.  3.3:  Outer  mean  velocity  on  the  body  sur¬ 
face  as  a  function  of  s. 


Inner  mean  Sow 


The  continuity  and  Navier-Stokes  equations  in  the  viscous  boundary  layer  in  optimal  coor¬ 
dinates  are  obtained  by  the  asymptotic  analysis  with  Rh  2>  1.  They  are 


8U  dV 


(3.12) 
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du  du  1  aub  2  _  i  dp  i  d2u 

+  +  Ub  U2  Rr  d^l  ’ 


(3.13) 


dP 

d ^ 


=  0, 


(3.14) 


Here  U  and  V  are  the  velocity  components  of  the  inner  mean  flow. 

Since  the  inner  mean  flow  is  two-dimensional,  a  stream  function  can  be  defined: 


where  p  is  the  inner  variable 


By  definition 


U  = 


™L  =  F' 

3TL 


(3.15) 


(3.16) 


(3.17) 


By  substituting  (3.17)  and  (3.18)  in  (3.13)  and  using  the  Bernoulli  equation 


(3.18) 


dP  _  ~  d  Ub 

d$L  bd$L’ 

it  is  obtained: 


F"'  +  FF"  +  2m($L)[l  -  F'2]  =  24>l(P'F4l  -  F*lF"),  (3.19) 

where 

dF  =  dF_  m  = 

dp'  *  34>l’ 

Equations  (3.12)-(3.14)  have  to  satisfy  the  following  boundary  conditions: 

U($L,  0)  =  0,  F($L,0)  =  0,  U (4>l,  oo)  =  1,  (3.20) 

which  are  the  no-slip  condition  at  the  wall  for  the  two  velocity  components  and  the  outer 
boundary  condition  which  states  that  the  velocity  must  be  equal  to  the  outer  velocity.  The 
boundary  conditions  for  equation  (3.19)  are  then: 

F(4>L,  0)  =  0,  F'($ L,  0)  =  0,  F'($l,  oo)  =  1,  (3.21) 

Equation  (3.19)  is  the  most  general  steady  mean  flow  equation,  in  optimal  coordinates 
and  the  flow  depends  on  <hL  therefore  it  is  not  similar.  It  is  numerically  solved  by  following 
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the  algorithm  in  ?  with  the  appropriate  boundary  and  initial  conditions.  In  section  4.1.3 
the  numerical  results  of  the  inner  mean  flow  are  discussed. 

Equation  (3.19)  is  parabolic  an  initial  condition  is  necessary.  The  initial  point  of  inte¬ 
gration  is  the  stagnation  point  where  the  outer  mean  flow  has  a  particular  form  in 
coordinates: 


Ub  = 


$L- 


Thus  in  this  particular  case  the  coefficient  m  in  equation  (3.19)  is  m  =  1/2  and  the 
equation  reduces  to  a  simpler  equation  with  similarity  because  the  function  F  does  not 
depend  on  <bL: 


F'"  +  FF”  -F'2  +  1  =  0,  (3.22) 

with  the  boundary  conditions 

F(  0)  =  0,  F'{  0)  =  0,  F'(oo)  =  1,  (3.23) 

Equation  (3.22)  is  the  well  known  Hiemenz  equation  which  is  numerically  solved  in  order 
to  have  the  behaviour  of  the  inner  mean  flow  close  to  the  stagnation  point. 

A  regular  perturbation  expansion  is  used  near  the  stagnation  point,  where  the  solution 
is  expanded  as  a  series  in  <f>L  and  terms  with  equal  powers  of  <hL  are  collected  . 


F  = 

Fo  +  (h,F,, 

(3.24) 

m(<3?  L 

)  =  \  +  C <hL. 

(3.25) 

The  leading  order  gives  the  Hiemenz  flow  equation  while  the  second  order  gives  an 
equation  for  F\  which  depends  on  Fq  and  it  is  numerically  solved. 

Fq"  +  FqFq  -  F'2  +  1  =  0,  (3.26) 

Ff  +  FqF"  +  F\Fq  -  2FqF(  +  2 C-  2Fq2  =  2 F^F[  -  2FI'F(/.  (3.27) 

The  total  mean  flow  is  given  by  the  composite  solution 

U*OT  =  F'*  +  U*  -U;.  (3.28) 

By  scaling  by  U£  the  following  is  obtained 

Utot  =  F  +  ^ - 1,  (3.29) 

U* 

and  dividing  and  multiplying  the  second  term  by  the  dimensionless  total  mean  flow 
is  given  by 
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,  u 

{Ttot  —  F  +  —  1. 

Ub 


(3.30) 


The  composite  solution  for  the  wall  normal  velocity  component  can  be  obtained  in  a 
similar  way. 


3.1.3  The  perturbation  flow 
Outer  perturbation  How 

The  upstream  vortical  flow  is  represented  by  a  convected  perturbation  of  a  uniform  mean 
flow  U^,  namely  the  flow  upstream  is 


V  =  i  +  eUoo^L  -tL,yL,zL),  (3.31) 

where  V  is  the  dimensionless  velocity  vector  with  components  u,  v  and  w  along  x,  y,  z 
and 


Uoo  =  6«>eifc(fciG*L+fc2GKL+fc3G*L-tL)  +  c.C-j  (3.32) 

with  k  =  2ttL*/\ kio  =  1,  k2c  =  AJ/A^,  k3c  =  AJ/A*  ,  AJ,A^,A*  are  the  wave 
lengths  in  the  streamwise,  normal  and  spanwise  directions  thus  kk\G  ,kk2G  ,kk3G  are  the 
dimensionless  wavenumber. 

The  solution  for  the  perturbation  flow  in  the  outer  region  is  of  the  form  suggested  by 
?,  where  the  flow  is  given  by  a  two-dimensional  steady  mean  flow  and  a  three-dimensional 
unsteady  small  amplitude  perturbation  flow: 


{u,v,w,p}  =  {^(xL,yiJ,F(a;L,;yL),0,P(zL)}  +  e{u,  v,  w,p}(xL,  yL)e*fc(fc3G*L  H  (3.33) 

The  solution  (3.33)  is  substituted  in  the  continuity  (3.1)  and  Navier-Stokes  equations 
(3.2)  to  obtain  the  following  linearised  Navier-Stokes  equations 


V  •  Ci  =  0, 


(3.34) 


— ^ —  +  U  •  Vu  +  u  •  VU  =  — Vp  +  —  V  •  (Vu),  (3.35) 

-Rl 

where  u  =  {u,  v,  w}  and  U  =  {U,  V}. 

The  general  solution  to  (3.35)  suggested  by  ?  is 


u  =  V(^L  +  uH,  (3.36) 

where  V(/?L  is  the  part  of  the  perturbation  velocity  related  with  the  pressure  fluctua¬ 
tions  and  uH  is  the  effect  of  the  imposed  upstream  vorticity  field.  This  is  called  Helmholtz 
decomposition  and  ?  proved  the  existence  and  uniqueness  of  this  solution  for  Navier-Stokes 
equations.  The  latter  part  of  the  perturbation  velocity  is  necessary  because  the  upstream 
boundary  condition  would  not  be  satisfied  only  by  the  potential  part  since  upstream  the 
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pressure  fluctuations  are  zero  thus  the  potential  is  zero  and  it  cannot  match  the  imposed  up¬ 
stream  perturbation  velocity.  The  velocity  uH  satisfies  the  homogeneous  equation  (without 
the  pressure  term).  In  the  case  of  a  flat  plate,  uH  coincides  with  Uqo  (initial  gust)  because 
there  is  not  distortion  due  to  the  curvature.  This  can  mathematically  be  observed  in  the 
equations  for  uH. 


<9uh 

dtL 


+  U  •  VuH  +  uH  •  VU 


^V-(Vuh). 


(3.37) 


The  above  equation  contains  the  mean  flow  U  which  for  a  flat  plate  is  uniform,  thus  uH 
for  a  flat  plate  can  only  be  equal  to  the  initial  gust,  used  as  upstream  boundary  condition. 
This  is  the  case  studied  in  ?  where  the  solution  in  region  I  is  given  by  (3.1)-(?).  The  velocity 
must  satisfies  the  continuity  equation  which  in  their  case  becomes  the  Laplace’s  equation 
(3.3)-(?)  for  the  perturbation  potential. 

For  the  case  studied  herein  the  mean  flow  is  given  by  a  uniform  flow  plus  a  source  flow 
and  therefore  its  gradient  is  not  zero.  Thus  there  is  a  distortion  in  the  perturbation  velocity 
held. 

From  the  continuity  equation  iph  is  given  by 


V2(/9l  =  —V  •  uH.  (3.38) 

The  pressure  p  is  given  by 

P  =  -  U  •  V(^L  +  (3-39) 

<76l  -TLl 

One  of  the  hypothesis  at  the  beginning  stated  that  A*  <C  L* .  Equation  (3.37)  was 
obtained  by  using  the  outer  scale  L*  and  1  / Rh  is  the  small  parameter.  WKBJ  theory  is 
used  to  obtain  the  approximated  solution  of  a  linear  differential  equation  whose  highest 
derivative  is  multiplied  by  a  small  parameter  (?).  In  this  case  the  small  parameter  used 
to  obtain  the  solution  is  absent  from  the  equations  but  the  scale  of  perturbation  motion  is 
much  smaller  than  the  characteristic  length  of  mean  how.  Moreover  the  spatial  scale  of  the 
perturbation  is  much  smaller  than  L*  which  is  used  in  analogy  to  the  application  of  WKBJ 
theory  to  classical  ordinary  differential  equations,  where  the  equations  are  initially  scaled 
by  the  larger  scale.  This  case  is  very  similar  to  that  one  of  ?  (page  223)  where  there  is  not 
a  small  parameter  in  the  equations,  then  it  is  not  obvious  how  to  use  WKBJ  theory.  By 
using  WKBJ  theory  the  solution  to  (3.37)  can  be  taken  of  the  form 

uH  =  eiks  (u0  +  ^ui  +  ...^  e*fc(fc3G*L-tL),  (3.40) 

where  k  =  2ttL*/X^  S>  1.  By  substituting  (3.40)  in  (3.37)  and  equating  coefficients  of 
inverse  powers  of  k,  the  following  are  obtained: 

1  —  U  •  Vs  =  —77— (Vs  •  Vs  +  fcfG),  (3.41) 

, _ .  2^  ^ 

U  •  Vuo  —  — — Vs  •  Vuo  —  — — uoV2s  +  Uo  •  VU  =  0,  (3.42) 
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where  R =  R^/k  =  0(1). 

It  has  to  be  remarked  that  in  obtaining  equation  (3.41)  the  following  assumptions  have 
been  made: 

•  AJ,A|,,A*  <  L*  namely  the  wavelengths  of  the  disturbance  in  the  three  spacial  direc¬ 
tions  are  all  of  the  same  order  and  they  are  all  smaller  than  the  characteristic  length 
of  the  body  L*. 

•  R\q  is  of  order  1.  This  assumption  permitted  to  use  WKBJ  theory  to  expand  the 
solution  for  the  velocity  as  in  (3.40).  This  theory  is  used  to  obtain  the  approximated 
solution  of  partial  differential  equations  whose  the  highest  derivative  is  multiplied  by  a 
small  parameter.  The  equation  for  uH  is  (3.37)  where  uH  is  a  function  of  only  xh  —  yL, 
the  V  operator  is  in  coordinates  xL  —  yL  and  fL  is  the  time  scaled  by  the  mean  flow 
velocity  and  L*.  The  small  parameter  in  this  equation  is  1  /Rh.  In  WKBJ  theory 
the  solution  expands  as  in  (3.40)  where  the  phase  is  multiplied  by  a  large  parameter, 
which  should  be  the  inverse  of  the  small  parameter  in  the  equation  (in  this  case  the 
Reynolds  number  Rh).  We  actually  choose  to  multiply  the  phase  by  k,  following  ?  . 

It  is  convenient  to  use  the  optimal  coordinates.  Here,  in  the  outer  region,  the  potential 
and  stream  function  used  are  $L  and  \kL  associated  with  the  outer  mean  flow  U.  A  new 
function  is  introduced 


sol  =  ~^—s  +  i(<J*L  -  4>l0), 


R 


A<£ 


where  <hLo  is  an  arbitrary  constant  and  sql  =  Sq/L*.  Then  (3.41)  becomes 


(3.43) 


(3.44) 


where 


q2  = 


4  i 


2  _l_  ik3 G 


R 


- 1, 


with  |U|  =  |V3>L|.  Equation  (3.44)  is  called  Eikonal  equation. 

The  boundary  conditions  on  all  the  four  boundaries  are  taken  as  equation  (3.27)  in  ? 
which  is  the  solution  for  equation  (3.41)  for  a  constant  mean  flow,  thus  it  can  be  used  as 
the  solution  far  from  the  body,  where  the  mean  flow  is  uniform.  The  solution  is 


s  =  k2GVc  + 


\nl,Q2  -  *!o  - 


(xL  -  XL0). 


(3.45) 


By  substituting  Q2  for  |U|  constant  in  the  above  equation  the  following  is  obtained: 


s  =  k2GVc  + 


iRu  klc  -  lRM'1 


(xL  -  xL0). 


The  problem  is  simplified  by  considering  an  inviscid  outer  perturbation  flow.  Under  this 
assumption  new  hypothesis  can  be  introduced: 
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1.  S*  <C  A*,  where  5*  is  he  boundary  layer  thickness. 

2.  A*  <  AJ  <  L*  and  A*/A|  =  <D(k*),  A %/L*  =  0{k *),  A*/L*  =  0(k%),  \*J\%  =  0(1) 
where  k$  is  the  dimensionless  streamwise  wavenumber  given  by  27rA*/AJ,  which  under 
this  hypothesis  is  small. 

Each  hypothesis  has  an  important  physical  meaning. 

•  S*  <C  A*.  This  ratio  between  the  boundary  layer  thickness  and  the  spanwise  wave¬ 
length  is  taken  as  in  region  II  of  ?.  This  first  implies  that  there  is  not  displacement, 
like  in  region  I  of  ?,  thus  the  streamlines  are  only  deflected  by  the  actual  body  and 
not  by  the  additional  displacement  due  to  boundary  layer.  If  displacement  effect  was 
considered  the  streamlines  external  to  the  boundary  layer  would  be  deflected  upward 
of  a  distance  which  is  exactly  the  boundary  layer  thickness.  The  displacement  effect 
gives  rise  to  the  concept  of  an  effective  body,  namely  the  free  stream  sees  the  contour 
of  the  boundary  layer  as  if  it  was  the  body.  In  this  region  the  outer  flow  is  still  inviscid, 
which  is  in  fact  the  case  that  is  considered  here  and  the  boundary  region  equations 
simplify  to  the  boundary  layer  equations. 

•  A*  <C  AJ.  This  ratio  between  the  two  wavelengths  in  the  spanwise  and  streamwise 
directions  is  the  disturbance  hypothesis,  also  taken  in  ?.  It  was  experimentally  demon¬ 
strated  that  the  disturbances  associated  with  laminar  streaks  evolves  over  distances 
which  are  larger  in  the  streamwise  wavelength  than  in  the  normal  and  spanwise  wave¬ 
lengths.  Under  this  assumption  the  dimensionless  streamwise  wavenumber  is  small 
and  during  the  analysis  of  the  inner  region  it  will  be  seen  that  this  allows  to  obtain 
parabolic  equations  inside  the  boundary  layer. 

•  A |<L*.  This  ratio  between  the  characteristic  length  of  the  body  and  the  streamwise 
wavelength  allowed  to  use  WKBJ  method  during  the  analysis  of  the  outer  flow,  as 
already  stated. 

After  all  these  considerations,  the  solution  to  the  Eikonal  equation  is  the  (3.26)  in  ?, 
which  is  valid  in  the  inviscid  limit 


sol 


2 

-Raj. 


k'2G  4'  L  + 


+  i($L 


^lo)- 


(3.46) 


It  can  be  easily  demonstrated  that  the  solution  3.46  satisfies  both  the  Eikonal  equation 
3.44  and  equation  3.41  in  the  limit  of  R\^  — >•  oo  and  in  a  region  far  from  the  stagnation 
point,  where  | U| 2  is  of  order  1.  In  fact  the  Eikonal  equation  becomes 


f  ds0L\2  (  Osql  A  2  _ 

and  sol  in  the  inviscid  limit  is 

■sol  =  ^(4>l  —  3>lo), 


(3.47) 
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which  satisfies  3.47.  Equation  3.41  in  the  inviscid  limit  of  R\ #  — >•  oo  and  in  $L  —  TL 
coordinates  becomes 


i  —  iur 


ds 


=  0, 


(3.48) 


because  the  inviscid  outer  mean  flow  is  parallel  to  the  streamlines  (its  only  component 
is  along  <h).  Equation  (3.48)  is  satisfied  by  the  solution 


s  =  fc2GTL  + 


1 


|U| 


dK- 


(3.49) 


The  streamwise  coordinate  to  be  used  is  $  =  which  is  the  streamwise  coordinate 
scaled  by  the  streamwise  disturbance  wavelength.  The  right  scale  for  the  disturbance  is  4> 
because  the  disturbance  associated  with  laminar  streaks  evolves  over  distances  which  are 
larger  in  the  wavelength  along  4>*  than  in  those  ones  along  T*,  z*.  The  4>*  coordinate  is 
therefore  scaled  by  AJ.  This  can  be  better  understood  by  remembering  that 

u-£, 

and  since  x  is  rescaled  by  AJ,  also  the  potential  function  must  be  rescaled  with  the 
same  quantity  in  order  to  have  the  velocity  of  order  1.  Time  is  scaled  in  the  same  way 
because  streaks  modulate  slowly  in  time.  The  normal  coordinate  T  is  scaled  by  the  shortest 
disturbance  wavelength  A*.  The  solution  for  the  phase  factor  in  the  new  scales  is 


s  =  + 


1 


4  lui 


(3.50) 


where  is  the  dimensionless  normal  wavenumber. 

Once  the  solution  for  sql  (then  for  s )  is  obtained,  a  solution  for  uo  can  be  found.  The 
coordinates  are  changed  as  for  the  phase  s  and  the  following  is  obtained: 


U  •  |U|Vu0  -  rJ^-|U|2Vs  •  Vu0  -  — u0|U|2V2s  +  u0  •  |U| VU  =  0,  (3.51) 

where  the  V  operator  is  in  (T^Tl)  coordinates  and  the  boundary  conditions  are 

u0eiks  — >  u°°eifc(fcl^L+fc2GyL)  as  TL  — »  oo,  (3.52) 

u0eiks  ->  u°°eik^1GXh+k2Gy^  as  TL  ^ -oo.  (3.53) 

WKBJ  theory  is  employed  again  to  find  the  solution  for  the  perturbation  potential: 

^=e~Y  ( <A)  +  (3.54) 

and  tpo  is  found  by  substituting  the  above  equation  into  the  continuity  equation  (3.38): 
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_i u°°  •  (Vs  +  k3Gk ) 

</?0  V2s  +  k\G 

where  k  is  the  unit  vector  along  the  spanwise  direction. 

The  total  outer  perturbation  velocity  is  found  by  equation  (3.36). 
The  outer  perturbation  pressure  is  given  by 


(3.55) 


P  =  pe(.TL,yL)elfc(fc3G2:L  *l),  (3.56) 

and  it  is  found  by  the  outer  perturbation  potential  as 

P  =  -  u  •  V<^L  +  v^vVl.  (3.57) 

-ti  l 

The  above  formula  is  shown  to  be  valid  by  ? .  The  outer  perturbation  pressure  is  of  order 
k$  as  the  inner  perturbation  pressure  (leib-wundrow-goldstein-1999).  This  is  an  interesting 
result  because  it  means  that  the  inner  pressure  matches  the  outer  pressure  .  The  pressure 
is  written  as 


pe  =  k<s>po- 

where  po  is  of  order  1  and  it  is  given  by 

p0  =  ^o-|U|2^g.  (3.58) 

09 

3.1.4  Inner  perturbation  flow 
Boundary  region  equations 

The  solution  in  the  inner  region  is  of  the  form  suggested  by  ?: 


{w<p,wy,wz,p}  =  {U($,r)),V($,r)),0,P($)}+e  {  u0, 


2 . 


ikz 


R 


-V0,W0 


’  \JRk<$, 


|  ($,??K(fcz; 

(3.59) 


where  d>  is  of  order  1,  r)  is  of  order  1,  kz  =  27t  is  the  spanwise  wavenumber,  U,  V 
are  the  components  of  the  2-dimensional  mean  flow,  uojUojleo  are  the  components  of  the 
3-dimensional  perturbation  flow.  Navier-Stokes  equations  in  optimal  coordinates  are 


~  dwo>  ~  divw  8wz 

U — -  +  U — -  H - - 

dT  9T  dz 


(3.60) 


1  dp  1 

iji  a$  +  r 


u 2 


dw$  dw$  1 
dt  d<f>  5T  U 


_  _  „„  1  d2wq>  5  In  U  f  duuv  dw$\  2  a  In  U  dwz 

a<f>2  a^2  jj2  dz2  a ^  \  a<i?  aT  )  u  a$  dz 


a2  re#  d2w<s> 

I  T  <1  I 


(3.61) 
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1  dp  1 
+  R 


1  dwq,  dwx ir  dwq,  1  dwqj  d  log  U  ,  9  9  . 

+  ”*^+--8f  +  “W  W  +  = 


u 2 

|2'>"-  <92xi;^ 


d2w^ 


+  fl^2  +  jj2  qz  2 


+ 


<9$ 

1  <92u><i> 


o  <9  In  U  f  duhj,  dw$\  2  d  In  U  dwz 

+  2  d<S>  \~d¥  _!)¥)-  u~d^  37  ’ 

(3.62) 


1  dp  1 
IJ3  dz  +  R 


d2 


Wz 


1  dwz  dwz  dw z  1  dwz  (  d  log  U  d  log  U 

[/2  dt  <9$  ail/  u  dz  y  <94>  <9' P 

d2wz  1  d2wz  wz  f d2JJ  d2u\  o  fdwzd\nU  dwzdlnU 

~ 1"  O.t.O  ~f"  O  O  d~  '  '  [  r\  t  O  ~f"  ^I.T.O  1  ~ 1"  "I  O  d~ 


<94> 2  dT2  jj2  Qz 2  jj  \  <94>2  dx V2 


<94>  <94» 


d T  <9T 
(3.63) 


By  substituting  (3.59)  in  (3.60)-(3.63),  taking  the  outer  inviscid  mean  velocity  along  the 
body  surface  and  putting  together  terms  of  order  e,  the  following  are  obtained: 


duo  _  duod_  dv o  w0  ikz  _ 

<94>  dp  2<&  dp  ub  k(j> 

Since  in  the  boundary  region  kz  is  of  order  1,  uq  and  Do  need  to  be  rescaled: 

ikz 


ikz  _ 
u0  =  —u, 

rC<|> 


_ 

=  UV' 

Av<|> 


Wo  =  w, 

where  u,  v  and  id  are  now  the  functions  to  find.  Thus  (3.64)  becomes 


(3.64) 


(3.65) 


du  du  p  dvw 
<94>  dp  24*  dp  jjb 


(3.66) 


iu  dF  du  _  d2F  _d2F  p  dF  du  F  du  _d2F 

+  TTTT+M^ - “TTWrV  W  W - ^  +  V  ~ 


1  **-£*.  (3.67) 


U2  dp  <94>  d<bdp  dp2  2<J>  3$  dp  2$  dp  dp2  2$  dp2  jj? 


iv  dF  dv  F  dv  1  <9  /  dF\ 
U 2  dp  <94>  24>  dp  24>  dp  \  dp  ) 


_  f  1  [  d  (  dF 

~u\  (2$)2  i'&p  \~dp 


-  F 


+ 


d2F 


dF  dv  _  d2F 
<9<E>  dp  1  d&dp 


11  dp 

ZZZc^  <.  i  ■  1  I 


d<s>2 
1  d2v 


p  d2F  1  dF 

4>  d&dp  4*  94> 


1  2 

—  ^ ~K  V, 


u2  24>  dp  24>  dp2  JJ? 


(3.68) 


iw  dF  dw  dF  dw  F  dw  k2  _  d2w  1  1  2  - 

u2  dp  <94>  <94>  dp  2$  dp  jjf  dp2  24-  [/2  '  v  '  ; 

In  (3.67),  (3.68)  and  (3.69)  terms  of  the  smallest  power  of  kq>  are  put  together  and 
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kz 

\/  Rk<p 


0(1). 


Equations  (3.66)-(3.69)  need  7  boundary  conditions  in  r/.  Three  of  them  are  the  no  slip 
conditions  for  u,  v  and  w,  namely 


u  =  v  =  w  =  0,  at  i]  =  0.  (3.70) 

The  streamwise  perturbation  velocity  u  inside  the  boundary  layer  is  of  order  l//c<j>  (see 
(3.59)),  while  the  streamwise  perturbation  velocity  u  outside  the  boundary  layer  is  of  order 
1  (see  (3.33)).  Since  u  is  larger  than  u  it  can  not  match  with  the  outer  solution  for  77  — >•  00, 
then 


u  — ^  0  as  — >•  00.  (3-71) 

The  perturbation  velocities  v  and  w  can  match  the  outer  perturbation  solutions  v  and 
w  because  they  are  both  of  order  1,  then 


v  — >•  —  — T-j  as  rj  — ^  00,  (3.72) 

kUW 2$ 

w 

w  — >  —  as  ii  — ^  00.  (3.73) 

Ub 

The  order  of  the  outer  perturbation  pressure  around  the  Rankine  body  is  the  same  of 
the  inner  perturbation  pressure,  thus  at  large  distance  from  the  body  the  inner  pressure 
matches  the  outer  perturbation  pressure: 


/  R  \  1/2  _ 

P  \  1 —  Po  as  r]  — >•  00. 

ik  \  k<\, ) 


(3.74) 


The  initial  conditions  for  boundary  region  equations  are  valid  in  the  region  where  $<1 
and  they  are  obtained  in  the  form  of  a  power  series. 
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4.1  Mean  flow  solution 

4.1.1  Mean  pressure 
The  mean  pressure  coefficient  is  defined  as 

Cp  =  1  -  U2. 

The  pressure  changes  along  the  body  surface  because  of  the  curvature  and  it  is  studied 
as  a  function  of  s.  The  behavior  of  Cp  is  shown  in  figure  4.1.  At  the  stagnation  point  Cp  =  1 
and  since  there  is  symmetry  the  slope  of  the  curve  is  zero.  Then  the  fluid  accelerates  and  the 
pressure  decreases  until  the  fluid  starts  to  decelerate  and  the  pressure  increases.  This  result 
is  also  obtained  by  ?  and  his  curve  and  the  one  numerically  obtained  perfectly  overlap  as 
shown  in  figure  4.1. 


Fig.  4.1:  Mean  pressure  gradient  as  a  function  of 
s  along  the  body. 

4.1.2  Boundary  layer  thickness 

The  boundary  layer  thickness  is  the  distance  from  the  body  to  a  point  where  the  velocity 
reaches  the  free  stream  velocity  value,  in  particular  it  is  defined  as  the  distance  where 
u  =  0.99C/oo-  This  distance  can  be  numerically  found  by  the  definition  of  r/,  which  is  the 
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inner  variable  used  to  find  the  solution  in  the  viscous  boundary  layer.  In  optimal  coordinates 
r]  is  define  as 


V  = 


Rn 

2$t 


Once  obtained  the  solution  of  the  inner  mean  flow  (next  section)  the  value  of  rj  where 
the  velocity  has  the  free  stream  value  can  be  found  for  every  streamwise  coordinate  $L  and 
this  value  is  supposed  to  be  different  on  the  curvature  where  the  solution  is  not  similar. 
Two  important  checks  can  be  done  analytically  and  numerically:  at  the  stagnation  point 
and  downstream,  where  the  flow  is  uniform,  the  boundary  layer  thickness  is  well  known.  In 
particular  at  the  stagnation  point  the  potential  and  stream  function  are: 


T  27T  ,  , 
=  -^LyL, 


$T,  = 


TTX^ 

A 


where  the  coordinates  with  prime  are  the  ones  used  at  the  stagnation  point  (figure  3.1). 
By  substituting  the  two  above  equations  in  the  definition  of  r/  the  following  is  obtained: 


V  =  Vl 


27 tRl 
A 


Thus  the  boundary  layer  thickness  for  the  stagnation  point  is  constant  and  equal  to: 


Vc  =  '7<5 


A 


2itRt 


where  rjs  is  the  value  of  i]  where  the  velocity  is  0.99  and  it  is  numerically  found  to  be 
r]s  =  2.4.  This  result  is  in  agreement  with  ?: 


where  a*  =  cdJ^/L*,  a  =  27r/A.  By  substituting  the  value  of  a*  in  the  boundary  layer 
thickness  found  by  ?  it  is  obtained: 


8t.  =  2.4 


A 


2vr Rr  ’ 


which  is  the  thickness  previously  found.  The  same  analytical  check  can  be  done  down¬ 
stream  where  the  potential  and  stream  fnction  are: 


=  Vl, 


$l  =  3A, 


and  therefore 


2/l  =  VS 


2xl 

R, 


1/2 

Downstream  the  boundary  layer  grows  as  xL'  .  On  the  body  is  no  possible  to  find 
analytically  the  boundary  layer  thickness  since  the  expressions  of  the  potential  and  stream 
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function  are  more  complicated.  Therefore,  the  boundary  layer  thickness  is  numerically 
evaluated  and  its  value  is  compared  with  Hiemenz  boundary  layer  at  the  stagnation  point 
and  Blasius  boundary  layer  downstream. 

The  behavior  of  the  boundary  layer  thickness  along  the  body  is  shown  in  figure  4.3.  At 
the  stagnation  point,  the  Hiemenz  boundary  layer  is  constant  because  there  is  a  perfect 
balance  between  the  wall-normal  viscous  diffusion,  which  lets  the  boundary  layer  thickness 
increase  and  the  favourable  pressure  gradient,  which  lets  the  boundary  layer  thickness  de¬ 
crease.  Downstream  the  pressure  gradient  effect  never  overcomes  the  wall-normal  diffusion 
and  eventually  wall-normal  diffusion  takes  over  completely. 

4.1.3  Inner  mean  velocity 

In  this  section  the  finite  difference  solution  of  the  steady  mean  flow  equation  in  optimal 
coordinates  is  presented.  Figure  4.2  shows  the  composite  solution  for  the  total  mean  flow 
given  by  equation  (3.30)  at  s  =  1.48.  In  figure  4.3  the  profiles  of  the  velocity  are  shown 
at  different  positions  on  the  body.  At  each  position  the  velocity  reaches  a  different  free 
stream  value,  as  expected  since  the  outer  mean  velocity  is  not  constant.  The  point  where 
the  velocity  has  the  free  stream  value  increases  following  the  boundary  layer  thickness. 


Fig.  4.2:  Composite  solution  of  the  mean  flow. 


4.2  Outer  How  solution 

Figure  4.4  shows  the  outer  perturbation  potential  at  different  streamlines.  It  is  the  con¬ 
tribution  to  the  total  perturbation  velocity  due  to  the  body  therefore  it  is  zero  at  large 
distance  from  the  body  (T  =  5).  Its  behaviour  changes  by  approaching  the  body.  Along 
the  streamline  T  =  1  the  potential  is  zero  upstream  and  downstream  where  the  mean  flow 
is  uniform  and  it  increases  on  the  body  where  the  mean  flow  accelerates.  This  behaviour 
is  more  evident  along  the  streamline  closer  to  the  body  (T  =  0.09),  where  the  potential 
increases,  then  it  rapidly  decreases  and  it  increases  again  to  reach  a  constant  value  which 
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Fig.  4.4:  Outer  perturbation  potential. 


i 

i 


i 

i 

i 


/ 


Fig.  4.5:  Outer  perturbation  pressure. 
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Fig.  4.6:  Outer  perturbation  velocity. 


is  kept  downstream  where  the  flow  is  parallel.  Figure  4.5  shows  the  outer  perturbation 
pressure  which  is  related  to  the  outer  perturbation  potential.  The  pressure  is  zero  at  large 
distance  from  the  body  ('k  =  5),  but  it  changes  by  approaching  the  body.  Along  the  stream¬ 
line  \k  =  1  it  is  zero  upstream  and  downstream  where  the  mean  flow  is  uniform  and  it  varies 
when  the  flow  accelerates.  At  \k  =  0.09  the  pressure  is  zero  upstream,  then  it  varies  on 
the  body  and  it  is  constant  downstream  where  the  flow  is  uniform.  Figure  4.6  shows  the 
behavior  of  the  spanwise  outer  perturbation  velocity  component  at  different  streamlines. 
Far  from  the  body  it  has  the  same  constant  value  as  the  initial  gust.  At  T  =  1  it  is  constant 
upstream  and  downstream  where  the  mean  flow  is  uniform  and  it  decreases  by  approaching 
the  body.  This  behaviour  is  more  evident  at  *k  =  0.09  closer  to  the  body.  The  other  two 
velocity  components  have  a  similar  behaviour. 
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4.3  Boundary  region  solution 

Figures  4.7  and  4.8  show  the  amplitude  of  the  streamwise  and  spanwise  velocity  components 
at  different  values  of  the  streamwise  coordinate  4>.  The  effect  of  the  pressure  gradient  is 
to  modify  the  streaks  which  do  not  have  a  single  peak  as  in  the  flat-plate  boundary  layer. 
The  shape  of  the  streamwise  velocity  amplitude  is  different  from  the  results  of  ?  and 
this  is  more  evident  at  small  4*.  The  streaks  are  initially  amplified  and  then  they  decay 
downstream.  The  behaviour  of  the  spanwise  velocity  amplitude  is  due  to  the  free-stream 
matching  requirements.  The  effect  of  k  on  the  streaks  at  $  =  0.5  is  shown  in  figure  4.9.  As 
k  decreases  the  streamwise  velocity  amplitude  increases. 


Fig.  4.7:  Profile  of  streamwise  perturbation  ve-  Fig.  4.8:  Profile  of  spanwise  perturbation  ve¬ 


locity  u  at  the  indicated  values  of  <f>. 


locity  w  at  the  indicated  values  of  $. 


Fig.  4.9:  Profile  of  streamwise  perturbation  ve¬ 
locity  u  at  the  indicated  values  of  k 
and  #  =  0.5. 
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5.  CONCLUSIONS 


The  downstream  evolution  of  a  perturbed  incompressible  flow  at  high  Reynolds  numbers 
around  a  Rankine  body  has  been  analyzed.  Far  from  the  body  the  flow  is  decomposed 
in  a  mean  flow  and  a  perturbation  flow.  The  first  one  is  given  by  a  uniform  flow  which 
encounters  a  source  flow.  The  solution  for  the  outer  perturbation  flow  is  found  by  using 
WKBJ  theory,  which  returns  the  phase  and  the  amplitude  of  the  perturbation  velocity  and 
consequently  the  perturbation  pressure.  The  viscous  effects  in  the  outer  perturbation  flow 
are  neglected. 

In  the  unsteady  boundary  layer  region  the  mean  flow  satisfies  the  steady  nonlinear  mean 
flow  equations  without  similarity  while  the  perturbation  flow  is  governed  by  the  linearized 
unsteady  boundary  region  equations  in  new  suitable  coordinates.  The  boundary  region 
equations  satisfy  appropriate  boundary  conditions  derived  from  matching  with  the  solution 
in  the  outer  region  and  no-slip  condition  at  the  wall.  Initial  conditions  for  boundary  layer 
equations  have  been  derived  in  terms  of  a  power  series  and  by  considering  that  at  the 
stagnation  point  the  inner  mean  flow  is  the  well  known  Hiemenz  flow. 

The  effect  of  the  mean  pressure  gradient  and  the  streamwise  wavenumber  on  the  stream- 
wise  velocity  profile  has  been  analyzed.  It  has  also  been  demonstrated  that  the  inner  per¬ 
turbation  pressure  matches  the  outer  perturbation  pressure  which  has  a  very  important  role 
in  the  boundary  layer  dynamics. 


5.1  Future  work 

The  rigorous  mathematical  approach  used  in  this  work  could  be  useful  to  extend  the  results 
to  more  complex  cases.  Some  ideas  have  already  been  developed. 

•  Considering  the  streamwise  wavelength  of  the  same  order  of  the  characteristic  length 
of  the  body  to  generalise  the  results  of  the  outer  perturbation  flow. 

•  Studying  the  effects  of  viscosity  on  the  outer  perturbation  flow  and  how  these  effects 
change  the  solution  inside  the  boundary  layer. 

•  Studying  the  evolution  of  order-one  disturbances  by  solving  the  nonlinear  boundary- 
region  equations. 

•  Using  nonlinear  streaks  as  base  flow  to  study  secondary  instability. 

•  Studying  the  perturbed  flow  around  a  wedge  body.  In  this  case  there  would  be  Falkner- 
Skan  boundary  layer  with  self  similar  equations  along  the  entire  body. 
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A.  DERIVATION  OF  CONTINUITY  AND  NAVIER-STOKES  EQUATIONS  IN 

OPTIMAL  COORDINATES 


Since  the  body  is  not  a  flat  plate  the  analysis  of  the  outer  and  inner  flow  is  made  in  a  new 
coordinate  system.  In  order  to  fully  understand  the  significance  of  this  choice,  it  is  necessary 
to  make  a  regression  to  explain  how  a  general  change  of  coordinates  is  made.  Switching 
from  x  —  y  orthogonal  coordinate  system  to  an  another  generic  orthogonal  coordinate  system 
$  —  T,  the  following  relation  is  valid 

dx2  +  d  y2  =  h%d&2  +  h%  d'L2, 

where  /i<j>  and  hq,  are  the  transformation  coefficients. 

It  can  be  easily  demonstrated  that 


(9$)2  _  (03)2 
,2  '  a  y  >  v  ox  > 

®  _  (  d<&\2  (  &®_\2  ’ 

\dx)  V  dy  )  \dx>  ^dy> 


(&*\2 

i  2  _  1  dx  ' 

fI>  —  (d}&\2(d®\2 

\dx>  \dy) 


(&$  \2 

2  dy  >_ _ 

(d$\2(dV \2 ’ 
'  dx  *  v  dy  ' 


(A. la) 


(A.lb) 


Now  the  general  transformation  coefficients  can  be  calculated.  A  special  orthogonal 
coordinate  system  is  chosen,  where  <f>  is  the  potential  function  and  ’L  is  the  stream 
function  of  the  outer  mean  flow  U  of  components  U  and  V.  By  definition  the  following 
relations  are  valid: 


~  _  3$  _  S'L 
dx  dy  ’ 

dy  dx 


(A. 2a) 
(A. 2b) 


By  calculating  the  transformation  coefficient  from  (A. la)  and  (A.lb)  the  result  is  that 

1  1 


—  h'xji  — 


U| 


u2  +  v2 


Then 


dx2  +  dy2  =  /i|d<L2  +  h\ d'L2  = 


d$2  +  d^2 
IUI2 


(A.3) 
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Fig.  A.l:  Outer  potential  velocity  in  x-y  plane  and  plane. 

The  velocity  U  is  a  vector  in  the  $  direction,  because  by  definition  a  potential  flow 
is  along  constant  streamlines.  Thus  it  has  two  components  in  the  x-y  plane,  but  just  one 
component  in  the  T-T  plane  along  <f>.  This  means  that  the  absolute  value  of  the  velocity  is 
the  same  value  as  the  T-component,  namely 


|U|  =  U*. 

Figure  A.f  shows  the  just  explained  relationship,  which  will  be  used  when  the  coordi¬ 
nates  of  the  equations  will  be  changed. 

The  reason  why  it  is  convenient  to  change  coordinates  and  in  particular  to  choose  the 
potential  and  stream  function  as  new  coordinates  is  that  in  this  new  system  the  body  can 
be  treated  as  a  flat  plate,  because  the  stream  lines  follow  the  shape  of  the  body  as  shown 
in  Figure  A. 2. 
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4- 


Fig.  A. 2:  Change  of  coordinates. 


Now  it  is  important  to  see  how  gradient,  divergence  and  curl  operations  change  when  the 
coordinate  system  changes.  The  function  a  is  a  scalar  and  A  is  a  vector  with  components 
(Aq>,Aq,,Az)  (these  components  are  in  the  new  directions);  The  new  coordinate  system  is 
'k,  z,  the  new  directions  are  e<j>,  eq>,ez  and  the  transformation  coefficients  are  /i$,  hq,,  hz. 
Then  the  following  relations  are  valid: 


V- A  = 


dAx  dAn 


dx 


+ 


dA7 


1 


dy  dz  hqdiq,hz 


d(hq,hzA<&)  d(h$hzAq,)  d{hq>hqi  Az 

t.  T 


<94* 


<9T 


dz 


~da  - da  f  da  „  1  da  „  1  da  „  1  da 

Va  =  ld~x+Jdy+kd~z=e^M  +  ^l^d^  +  ezVzWz' 


V  x  A  =  i 


1 

h^hq/hz 


dAz  dA, 


dy 


dz 


dAz  dAz 


dx  dz 


+  k 


dAv  dAz 


dx 


dy 


(A. 4a) 
(A. 4b) 

(A. 4c) 


eq>hq>  ( d(hzAz^  _  \  ( d{hzAz )  _  d(Ji^A^)\  +  ~  h  f  d(hq,Aq,) 


<9T 


dz  J 


<94* 


dz  J 


V 


Relations  (A. 4)  are  used  to  obtain  continuity  and  Navier-Stokes  equations  in  the  new 
coordinate  system.  In  the  case  considered  here 

hq>  =  hq/  =  , 

IUI 


hz  =  1. 

In  order  to  derive  continuity  and  Navier-Stokes  equations  in  optimal  coordinates  namely 
from  ( x,y,z )  to  (<h,T,2:),  the  following  equation  is  used  to  write  the  term  V •  (VV)  in  another 
way: 


v  •  VV  =  V(V  •  V)  -  V  x  V  x  V. 

By  the  continuity  equation  the  term  V  •  V  is  zero.  Thus  (A. 5)  becomes 


(A.5) 


<9(/i$  A$) 
<94* 
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v  •  VV  =  -V  x  V  x  V. 

Moreover,  it  is  easy  to  demonstrate  the  following  relation: 

IVI2' 


V  •  VV  =  V 


-  V  x  V  x  V. 


(A.6) 


(A.7) 


By  substituting  (A.6)  and  (A.7)  in  the  continuity  and  the  following  Navier-Stokes  equa¬ 
tions 


V- V  =  0,  (A. 8) 

8V  1 

_  +  V.VV  =  -Vp+-V-(VV),  (A-9) 

where  R  is  Reynolds  number,  it  is  obtained 

-  V  x  V  x  V  =  — Vp  —  -)-V  x  V  x  V,  (A. 10) 

R 

which  is  another  form  of  the  Navier-Stokes  equation. 

By  equations  (A. 4)  it  is  easy  to  see  the  the  continuity  equation  in  optimal  coordinates 


~  dw$  ~  duhi,  8wz 

U — -  +  U — -  H - - 

9$  d dz 


(A.ll) 


where  U  =  |U|,  =  u$/U,  =  uy/U  and  wz  =  w/U,  with  u$  and  the  velocity 

components  in  the  4>  and  vk  directions  respectively. 

In  order  to  find  the  Navier-Stokes  equations  in  optimal  coordinates  the  following  passages 
are  made: 
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}} 


and  thus  equation  (A.  10)  in  scalar  form  is 
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where  u,v,w  are  the  velocity  components  in  the  directions  (not  in  x,y,z  directions). 

Rearranging  and  simplifying  it  is  obtained 
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B.  NUMERICAL  PROCEDURE  FOR  BOUNDARY  LAYER  EQUATIONS 


In  order  to  solve  the  boundary  region  equations  a  marching  procedure  along  3*  is  used.  The 
method  is  based  on  second-order  central  differences  in  r]  and  backward  differences  in  3*.  The 
grid  used  to  compute  the  pressure  is  staggered  compared  to  the  grid  used  for  the  velocity 
components  as  shown  in  Figure  B.l. 


d<3> 


O - o- 

bb  b 


OJ  +  1 

o  i*  +  1 
-•  i 
o  F 


h 


oj-1 


I  ^ 

Fig.  B.l:  Grid  used  to  compute  the  velocity  components  and  pressure. 


With  the  second-order  central  differences  method  a  function  /  and  its  first  and  second 
derivatives  are  approximated  as  follows: 


/  =  fi, 


(B.l) 


dj_  =  fj+i  -  fj- 1 

drj  2  h 


(B.2) 


d2/  =  fj+ 1  -  2 fj  +  fj_i 

dr )2  hi2 


(B.3) 


With  the  second-order  backward  differences  method  the  first  derivative  of  a  function  / 
is  given  by: 


<9/  _  afj  +  bfbj  +  cfbbj 
<93*  ~  d<F 


where  a,  b  and  c  are  assigned  values. 


(B.4) 
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The  grid  of  pressure  is  staggered,  thus  the  value  of  the  pressure  in  (3.69)  and  its  deriva¬ 
tive  with  respect  to  r]  in  (3.68)  have  been  expressed  in  terms  of  the  staggered  index,  namely 


=  Pj*+ 1  +Pj * 
2 


(B.5) 


dp  =  pj*+ 1  -pj* 

Op  h 

as  can  be  easily  seen  in  Figure  B.l. 

The  index  j  goes  from  0  to  iV  —  1,  thus  the  boundary  conditions  are: 


(B.6) 


Uj=o  —  Vj=o  —  Wj=Q  —  0, 

(B.7) 

Uj=N~  1  =  0, 

(B.8) 

.  u(d>,  N  —  1) 

:7V- 1  =  l  ~  -  r^, 

kUb(®,N-  l)v/2$ 

(B.9) 

w(<P,  N  —  1) 

Wj=N-l  =  ~  ,  j 

(B.10) 

Ub($,N-  1) 

Pj=N-l  =  Po{$,N  -  1.) 

(B.ll) 

The  linearised  parabolic  system  have  a  block  tridiagonal  structure  which  can  be  solved 
with  Keller’s  Box  method  .This  block  elimination  method  consists  of  two  sweeps,  one  forward 
and  one  backward.  At  each  step  j  the  unknown  functions  are  computed  from  recursion 
formulas. 
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